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The problem of anomalous scaling in magnetohydrodynamics turbulence is considered within the 
. . . framework of the kinematic approximation, in the presence of a large-scale background magnetic 

' field. The velocity field is Gaussian, cS-correlated in time, and scales with a positive exponent ^. 

, Explicit inertial-range expressions for the magnetic correlation functions are obtained; they are rep- 

resented by superpositions of power laws with non-universal amplitudes and universal (independent 
of the anisotropy and forcing) anomalous exponents. The complete set of anomalous exponents for 
^ , the pair correlation function is found non-perturbatively, in any space dimension d, using the zero- 

mode technique. For higher-order correlation functions, the anomalous exponents are calculated to 
O(^) using the renormalization group. The exponents exhibit a hierarchy related to the degree of 
0^ ■ anisotropy; the leading contributions to the even correlation functions are given by the exponents 

from the isotropic shell, in agreement with the idea of restored small-scale isotropy. Conversely, the 
small-scale anisotropy reveals itself in the odd correlation functions : the skewness factor is slowly 
decreasing going down to small scales and higher odd dimensionless ratios (hyperskewness etc.) 
dramatically increase, thus diverging in the r ^ limit. 
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I. INTRODUCTION 



I In cosmical objects, small-scale evolution of the magnetic field B often takes place in the presence of a strong 
. large-scale magnetic field B° . It is, for example, what happens in the solar corona where, in spite of the typical value 
of the sun magnetic field (« 1 Gauss) , fields as intense as « 500 Gauss can be observed in solar flares. These highly 
, energetic and large-scale events coexist with small-scale turbulent activity, finally responsible of the dissipation of 

■ magnetic field energy. Modelling the way through which energy is stored and then dissipated is, consequently, not an 
easy task. 

In Ref. 01, the following description is proposed: a large-scale axial, e.g., directed parallel to some vector z, magnetic 
field B° is assumed to dominate the dynamics in the z direction, while the activity in the transverse plane can be 
satisfactorily described as quasi-bidimensional. This picture allows reliable numerical simulations in two dimensions, 
^ ■ from which it appears clear that the magnetic field tends to organize in rare large-scale structures separated by 
^ narrow current sheets. Deep investigation of small-scale intermittency properties is still not permitted by lack of 
spatial resolution. 

An interesting question raised by this problem, beside structure formation, is related to the role played by large-scale 

■ anisotropy on the small-scale statistics. Indeed, this is quite a typical situation in turbulence, where almost every 
large-scale forcing is not isotropic. Here, instead of taking the restoration of local small-scale isotropy for granted, as 
in the Kolmogorov theory of turbulence (^^] , we analyze in detail the effects of anisotropic large-scale contributions 
on the small-scale magnetic fluctuations. 

A wide interest has been recently devoted to this issue P-p^. From the viewpoints of theoretical and numerical 
analysis, focusing on a small number of indicators some arguments are given in favor of the small-scale isotropy 
restoration in the Navier-Stokes (NS) turbulence ||ll|,|l^ • On the other hand, investigating a larger class of anisotropic 
indicators, footprints of small-scale anisotropy become manifest ||^. The scenario thus appears extremely faceted and 
needs of further investigations. 

Recently, clear evidences of persistent small scale anisotropy have been found in Ref. where the statistical 
properties of a scalar field advected by the non-intermittent NS flow generated in a two-dimensional inverse cascade 
regime are investigated. 

Two main goals motivate this paper. On one hand, we give details of the results presented in the Rapid Communi- 
cation p^ ] where the effects of anisotropy on scaling exponents of the two-point magnetic field correlations have been 
addressed in the framework of the kinematic magnetohydrodynamics (MHD) problem. Non-perturbative expressions 
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for the scaling exponents were derived and their universahty proved. Specifically, there arises a picture of a non trivial 
statistical behavior, where anisotropic fluctuations are organized in a hierarchical order according to their degree of 
anisotropy. Contributions belonging to shells of higher anisotropic index decay faster, and the isotropic contribution 
finally dominates. 

However, the dominance of the isotropic contribution in the scaling exponents does not imply that large-scale 
anisotropy is irrelevant for the small-scale magnetic statistics. A deep investigation focused on a larger number of 
statistical indices (that is focused on the proper anisotropy indicators) has to be performed in order to highlight 
the way (if any) through which large-scale anisotropy manifests itself at small scales. This is the second aim of the 
present paper. Specifically, in addition to the non-perturbative results for the two-point correlations, we present new 
results dealing with higher-order magnetic correlation functions. Being more specific, we exploit the field theoretic 
renormalization group (RG) to obtain the anomalous exponents for higher-order magnetic correlation functions at the 
first order in ^, the exponent entering into the velocity covariance. In particular, we evaluate the odd-order correlation 
function exponents, from which dimensionless ratios like skewness and hyper-skewness are calculated. As a result, 
in three dimensions, the former behaves at the dissipative scale as Pe^^/^" while the latter as Pe^^^^^, Pe being 
the Peclet number (i.e., the equivalent of the Reynolds number for the NS turbulence). Notice the opposite signs 
appearing in the scaling exponents. They are the signature of persistent small-scale activities. Indeed, the first index 
is weakly scale dependent while the second is even divergent at small scales (i.e., Pe — s- oo). Let us remark that to 
restore isotropy at small scales all such indices should decay to zero as Pe grows. 

The same general picture is found numerically in Ref. |16|| in the framework of the passive scalar advection by NS 
flows. In addition, the results are in qualitative agreement with the first-order analytic expressions for the anomalous 
exponents obtained in for the passive scalar advected by a synthetic velocity field. 

The paper is organized as follows. In Sec. ||, we give the detailed definition of the kinematic MHD Kasantzev- 
Kraichnan model, which describes the passive advection of the magnetic field by the Gaussian white-in-time, self- 
similar velocity field. In Sec. [II, the field theoretic formulation of the model is presented. It allows for the derivation 
of the closed exact equations for the response function and equal-time pair correlation function of the magnetic 
field. From the homogeneous solutions (zero modes) of the pair correlation equation, scaling exponents of the pair 
correlation function are determined. In Sec. pv| , these exponents are found non-perturbatively, for any ^ and space 
dimensionality d. In Sec. we discuss the UV renormalization of the model and derive the corresponding beta 
functions and RG equations. The latter possess an infrared (IR) stable fixed point, which establishes the existence 
of anomalous scaling for all the higher-order correlation functions. The inertial-range behavior of these functions is 
determined by the scaling dimensions of cert ain t ensor composite operators; they are calculated in Sec. VI to the first 
order in ^ (one- loop approximation). In Sec. VIl, we employ the operator product expansion to give explicit inertial 
range expressions for various higher-order correlation functions. The results obtained are reviewed in Sec. VIII, where 
a brief comparison with the passive scalar problem is also given. 



II. DEFINITION OF THE KINEMATIC MHD KASANTZEV-KRAICHNAN MODEL 

In the presence of a mean component B° (actually supposed to be varying on a very large scale ~ L, the largest 
one in our problem) the kinematic MHD equations describing the evolution of the fiuctuating part B = B{x) of the 
magnetic field are : 

dtBo,+v ■ dBa^ B ■ dva+ B° ■ dvo, + Kod^Bc,, a = l,---,d- (2.1) 

Here and below x = {t, x}, d = \da — d/dxa}, = dada is the Laplace operator, d is the dimensionality of the x 
space and v = v{x) is the velocity field. Both v and B are divergence free (solenoidal) vector fields: daVa = daBa = 0. 



Equation (2.1) follows from the simplest form of Ohm's law for conductive moving medium, j = a {E + v x B/c), and 
the Maxwell equations neglecting the displacement current: dtB/c+ dxE = Q, dxB^ Anj /c = and d ■ B — 0. 
Here c is the speed of light, j is the density of the electric current, a is the conductivity, and kq = /Ana is the 
magnetic diffusivity. 



The term B° ■ dva in (2.1) effectively plays the same role as an external forcing driving the system and being also a 
source of anisotropy for the magnetic field statistics. 

In the real problem, v obeys the NS equation with the additional Lorentz force term cx (c? x B) x B, which describes 
the effects of the magnetic field on the velocity field. The framework of our analysis is the kinematic MHD problem, 
where the reaction of the magnetic field B on the velocity field v is neglected. We assume that at the initial stages 
B is weak and does not affect the motions of the conducting fluid: it becomes then a natural assumption to consider 
the dynamics linear in the magnetic field strength [p7|. It is also noteworthy that in more realistic models of the 
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MHD turbulence the magnetic field indeed behaves as a passive vector in the so-called kinetic fixed point of the RG 
equations (see Refs. ||l^ , |l9|] ). 

For general velocity fields the well-known closure problem arises even for the kinematic model. This means that 
the equations of evolution for the single-time multiple-space moments such as {Ba{t, ri) ■ ■ ■ B\{t, r„)) are not closed. 
The situation changes for white- in-time random velocity fields. The physical choice of a real turbulent flow governed 
by the NS equation is then replaced by an incompressible, self-similar advecting field, with Gaussian statistics and 
rapidly changing (J-correlated) in time. This last property allows us to write closed equations for the moments of the 
magnetic field B and to perform analytical (both perturbative and non-perturbative) approaches to the d-dimensional 
problem. Indeed, in the presence of a white-in-time random velocity field, the solution is a Markov process in the time 
variable and closed moment equations, sometimes called "Hopf equations" , can be obtained in analogy to the passive 
scalar case ||2^. Such models have attracted enormous attention recently (see, e.g., Refs. pX|pl| -p5[| and references 
therein) because of the insight they offer into the origin of intermittency and anomalous scaling in fully developed 
turbulence. We also note that the isotropic version of the kinematic rapid-change magnetic model dates back to 1967 



(see Ref. M) and was studied by Refs. |27-30| 



More precisely, we shall consider a simplified model in which v{x) is a Gaussian random field, homogeneous, isotropic 
and white-in-time, with zero mean and covariance 

{v^ix)vp{x')) = S{t - <') X„^(r) (2.2a) 

with 

Ka^ir) = Do cxp[ik . r], r^x x' , (2.2b) 

where Pct^(k) = 5ap — kakp/k"^ is the transverse projector, k is the momentum, k = |k|, Dq > is an amp litude 



factor, and < ^ < 2 is a free parameter. The IR regularization is provided by the cut-off in the integral (2.2) from 
below at fc ~ m, where m = 1/L is the reciprocal of the integral turbulence scale; the precise form of the cut-off is 
not essential. For < ^ < 2, the difference 

Sc.fi{r) = Ka.f,{()) ~ Ka.p{r) (2.3) 

has a finite limit for m ^ 0: 

5„Mr)=i?r« [(d + e-l)<^o/3-e^] , (2.4) 

with 

-DoT{-i/2) 



D 



(4^)'^2«(d + Or(d/2 + C/2)' 



where r(- • •) is the Euler gamma function (note that D > Q). It follows from Eq. ( p.4| ) that ^ can be viewed as a kind 
of Holder exponent, which measures the roughness of the velocity field. In the RG approach, the exponent ^ plays 
the same role as the parameter e = A — d does in the RG theory of critical phenomena [|3l|| . The relations 

50 = i?o/«o = A« (2.5) 

define the coupling constant 50 (i.e., the expansion parameter in the ordinary perturbation theory) and the charac- 
teristic ultraviolet (UV) momentum scale A. 



III. FIELD THEORETIC FORMULATION OF THE MODEL. DYSON EQUATIONS FOR THE PAIR 

CORRELATION FUNCTIONS 



The stochastic problem (2T), ( p.2| ) is equivalent to the field theoretic model of the set of three fields $ = {B' , B, v} 
with action functional 



S'($) = B' -dtB - {v ■ d)B + {B ■ d)v + {B° ■ d)v + n^d^B - vK-^v/2. 



(3.1) 



The first five terms represent the Martin-Siggia-Rose action (see, e.g., Refs. |31-33|) for the stochastic problem J^) 
at fixed V, and the last term represents the Gaussian averaging over v; K^^ is the inverse integral operation for (|2.2b) 
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and B' is a solenoidal response vector field. In (3.1) and analogous formulas below, the required integrations over 
{i, x\ and summations over the vector indices are implied, for example, 

B'dtB = jdtj dxB'^{x)dtBa{x), vK^^v = JdtJdxJ dx' Va{t,x)K-^{x - x')vp{t,x'). 



The formulation (3J_) means that statistical averages of random quanti ties in the stochastic problem (2J), (2.2) 
coincide with functional averages with the weight expS'($). The model (B.l) corresponds to a standard Feynman 



diagrammatic technique with the triple vertex B' [—{v ■ d)B + {B ■ d)v\ ~ B'^B fjVjVapj with vertex factor 



7' 



(3.2) 



where k, p and q are the momenta flowing into the vertex via the fields B' , B and v, respectively. Strictly speaking, 
the vertex Vapj has to be contracted with three transverse projectors, but we omitted them in order to simplify the 
notation. In most cases, transversality of Vap-y with respect to all its indices will be restored automatically owing to 
the contraction with bare propagators. The latter in the frequency-momentum (w, k) representation have the form: 



■Pa/3(k), 



{B^{u;,k)B'0{-iu,~k))o = {B'^{cj,k)Bp{~u;,~k))* - 

(B„(w,k)B^(-c^,-k))o = (B°-k)2(B,(cc;,k)B;,,(-cc;,-k))o(«a'(cc;,k)«^,(-c^,-k))o(B;j,(a;,k)B^(-c^,-k))o, 
(B„(w, k)vi3i-cj, -k))o = (B° • k){B^{ij, k)B'^,i-ij, -k))o(t;a'(w, k)vi3i-u;, -k))o, 

{B'Ju;,k)B'f,i-ij,-k))o^O, (3.3) 



(2.2) 



and the bare propagator {vaV/j)Q is given by Eqs. 

The magnitude B° = \B°\ can be eliminated from the action (3.1) by rescaling of the fields: B — > B°B, B' 
B' /B° . Therefore, any total or connected Green function of the form (B{xi) ■ ■ ■ -B(x„) B'{yi) ■ ■ ■ B'{yp)) contains 
the factor of (5°)""^'. The parameter B° appears in the bare propagators (3.3) only in the numerators. It then 
follows that the Green functions with n — p < vanish identically. On the contrary, the 1-irreducible function 
{B{xi) ■ ■ ■ B{xn) B'{yi) ■ ■ ■ S'(yp))i_ii. contains a factor of {B°Y~'^ and therefore vanishes for n — p > 0; this fact 
will be relevant in the analysis of the renormalizability of the model (see Sec. 0). 

The white-in-time character of v permits to exploit the Gaussian integration by parts (a comprehensive description 
of this techniques can be found, e.g., in Ref. Q) to obtain closed, exact equations for the equal-time correlation 
functions of the field B. This strategy has been used in Ref. jl5|. Below we give an alternative derivation of the 
equation for the pair correlation functions based on the field theoretical formulation of the problem (see also Ref. |2^] 
for the scalar case). 

The pair correlation functions ($$) of the multicomponent field <& satisfy standard Dyson equation, which in the 
component notation reduces to the system of two nontrivial equations for the exact correlation function CQ/3(ti',k) 
{Ba{Lo, k)Bfi{—uj, —k)) and the exact response function Gap{w, k) — {BaiuJ, k)_B^(— w, — k)). The latter is independent 
of B° (see above) and thus can be written as Gap{uJ, k) — Pai3{k)G{LU, k) with certain isotropic scalar function G(ijJ, k). 
In our model these equations, usually referred to as the Dyson-Wyld equations (see, e.g., Ref. |^), have the form 



G-\u;,k)P^p{k) = [-iLo + Kofc2]Po/3(k)-Sf/(^,k) 



(3.4a) 



C„0(c^,k) = |G(L^,fc)|^ (B°-k)^(i;„(c^,k)z;0(-a>,-k))o-hSf/ (co,k) 



(3.4b) 



where {BaBp)^ is given in Eq. ( ^.3| ), ^ and ^ are self-energy operators represented by the cor respo nding 
1-irreducible diagrams; the other functions S** vanish identically. It is also convenient to contract Eq. (3.4a) with 
the projector Paj^ik) in order to obtain the scalar equation: 



G "'"(oj, k) = —iijj + Kpfc^ 



(3.5a) 



where we have written 



(^,fc)^sf;^(c.,k)F„Mk)/(d-i). 



(3.5b) 



The feature characteristic of the rapid-change models like (^J) is that all the skeleton multiloop diagrams entering 
into the self-energy operators ^ and ^ contain effectively closed circuits of retarded propagators {BB') and 
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therefore vanish; it is also crucial here that the propagator (vv) in Eq. ( ^.2a ) is proportional to the S function in 
time. Therefore the self-energy operators in are given by the one-loop approximation exactly and have the form 



(3.6a) 



^B'B' 



(3.6b) 



The solid lines in the diagrams denote the exact propagators {BB') and {BB)\ the ends with a slash corre spon d to 
the field B', and the ends without a slash correspond to B\ the dashed lines denote the velocity propagator (2.2); the 
vertices correspond to the factor (3.2). The analytic expressions for the diagrams in Eq. ([3.6|) have the form 



Zj [lu, k) — ' 



(d-l) 



k) f did' f dq / \ / / X 



( k, p, q) 



(3.7a) 



Sf0^'(^,k) 



duj' 

'2tt j {2tt 



d'i n \ I I \ -^0 ^aia2(q) xr / , ^ 

— J V^a3ai(k,p,q)Ca3a4(w ,p) ^^^^.^^(-k, -p,-q), 



(3.7b) 



where k -I- q + p = 0, the vertex Vc/j^ is defined in Eq. (3.2), and the explicit form (2.2) of the velocity covariance is 
used. We also recall that the integrations over q should be cut off from below at q = m. 

The integrations over w' in the right hand sides of Eqs. (3.7) give the equal-time response function G{q) = 
(l/27r) J dui' G{uj' ,q) and the equal-time pair correlation function Cap{ci) — (l/27r) J dui' Ca0{uj' ,q); note that both 
the self-energy operators are in fact independent of lu. The only contribution to G(g) comes from the bare propagator 
( |3.3[ ), which in the t representation is discontinuous at coincident times. Since the correlation function (2.2a), which 
enters into the one-loop diagram for ^ , is symmetric in t and t' , the response function must be defined at t — t' 
by half the sum of the limits. This is equivalent to the convention 



G(q) = (l/27r) / du' {-iuj' + Koq^) 



,2\-l 



1/2 



and gives 



V-B'B/ ,N -Pa/3(k) f dq \ u I- N ^0-Paia2(q) Tr / , > 

^ ^) = —1 I 77;ZT7^°"3ai (k, p, q) Pa,04 (P) —i ^0.41302 (^k, -p, -q). 



2{d- 1) J {2TrY 



Substituting Eq. ( |3.2D into Eq. (3.S) after lengthy but straightforward calculation gives 



{Lo,k) = {-l/2)k^kpDo 



{2n)'i q'i+i 



The integration over q in Eq. (3.9) is performed explicitly using the relation 



dqfiq)^ 
q^ 



-f I dqfiq) 



and gives 



where we have written 



:^'^i.,k) = -k^^^^^pljim), 



J{m) 



dq 1 

{2t:Y '^f+i 



(3.8) 



(3.9) 



(3.10) 



(3.11a) 



(3.11b) 



Here and below Q = Sd/{2T:Y and 5^ = 27r'^/Vr(d/2) is the surface area of the unit sphere in d-dimensional space; 
the parameter m has arisen from the lower limit in the integral over q. 
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Equations (3.5), (3.11) give an explicit exact expression for the response function in our model; it will be used in Sec. 
for th e ex act calculation of the RG functions. Like in the scalar case, the exact response function differs from its bare 
analog (3^) simply by t he su bstitution kq ^ kq + ^Dq (d — 1) J{m)/2d. Below we use the intermediate expression (^). 
The integration of Eq. ( 3.4b ) over the frequency oj gives a closed equation for the equal-time correlation function; it 
is im por tant here that the lo dependence of the right hand side is contained only in the prefactor |G(w, Using 
Eq. (3^) the equation for Capi^) can be written in the form 



2(Kofc' + ^)C„;3(k) = (B° • k)2(i,„(c^, k)i;0(-a>, -k))o 



+ 



(2^ 



(k,p,q)C. 



(-k, 



(3.12) 



and using Eqs. (3_^) and (3^) it can be rewritten as 



dq Do 



{2-kY q'^+^ {'7aigQ2Caia2(p)-Pa/3(q) - Pai ^Qg Caa2 (P)^ai /3 (q) " Pca ^ai Cai/3 (p)PQa2 (q) } 
j2nY |PaiPa2'-'a/3lPj " ^ Oq/3 (^Kj j- . 



(3.13) 



For < ^ < 2, Eq. ( |3.13| ) allows for the limit m — > 0: the first three integrals in its right hand side are separately 
finite for m — 0; the last inte gral is finite owing to the subtraction, which has come from the contribution with ^ 
in the left hand side of Eq. (3.12). Indeed, the possible IR divergence of this integral at q = is suppressed by the 
vanishing of the expression in the curly brackets. In what follows we set m = 0. 

Equation (3.13) can also be rewritten as a partial differential equation for the pair correlation function in the 
coordinate representation, Capir) = {Ba{t, x)B/s{t, x + r)) [ we use the same notation Cap for the coordinate function 
and its Fourier transform]. Noting that the integrals in Eq. ( 3.13) involve convolutions of the fu nctio ns Cap(k) and 
DoPap{k)/k'^~^^, the Fourier transform of the spatial part ( 2.2b ) of the velocity correlation function ( |2.2| ), and replacing 
the momenta by the corresponding derivatives, ipa — > da and so on, we obtain: 



idaiSaa2) i.doi2^ctif3) Saia2 ^0.1^0.2 ^af3- 



(3.14) 



Note that the correlation function (2.2) enters into Eq. (3.14) only through the function Saia2 from Eq. (2.4), or, 
in other words, through the difference (2.3), which has a finite limit at m = 0. The m dependent constant part of 
( ^.2b ) vanishes under the differentiation in the first four terms in the right hand side of Eq. (3.14), and in the last 
term it is subtracted explicitly, owing to the subtraction in Eq. ( 3.13 ). Equation ( 3.14 ) should be augmented by the 
solenoidality condition: 



daCaB — 0. 



(3.15) 



For the nonstationary state, the function Cgg jt, r) = (Ba(t, x)Bg{t,x + r)) depends explicitly on t, and the term 
dt Cap appears on the right hand side of Eq. (3.14), see, e.g., Ref. |15||. 



IV. NONPERTURBATIVE RESULTS FOR THE SCALING EXPONENTS OF THE 2-POINT 

MAGNETIC CORRELATION FUNCTION 

In this section we focus our attention on the inertial-range behavior of the second-order equal-time correlation 
function Cap(t,r) = {Ba(t,r)Bp{t,0)) in the statistically steady state. As shown in Ref. jl^, a steady state is 
present when ^ < 1, ^ = 1 being the threshold of instability. As such threshold coincides with that of the isotropic 
problem [ p7[ , it follows that dynamo effect is thus not swithched on by anisotropic contributions. 

In the isotropic case, the analytic expression for the scaling exponent of Cap has been obtained in Ref. . It was 
also shown by the author of [ p7| that the anomalous exponent is universal, and the anomaly is associated with zero- 
mode solutions of the equations satisfyed by Cap- Higher-order correlation function exponents have been calculated 
to 0(0 in Ref. H by exploiting the RG. 

With respect to Ref. j27|, the main technical difference is that, in order to extract the anisotropic contributions to 
the isotropic scaling, the angular structure of zero modes has now to be explicitly taken into account. 
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To start our analysis, let us consider the closed Eq. ( [3.14 ) for Cap- For what follows, it is worth emphasizing two 
properties of Ca/B- 

(i) because of homogeneity, Caf3 is left invariant under the following set of transformations: 



r 



and 



(4.1) 



(ii) Caf3{r) = Capi—r), as it follows from (3.14) after the substitution r ^ — r. 

In the presence of anisotropy, the most general expression for the two-point magnetic correlations, Capir), in the 
stationary state involves five (two in the isotropic case) functions depending on both r = \x — x'\ and z = cos 6* = 
B ■ r/r, where B is the unit vector corresponding to the direction selected by the mean magnetic field. Remark 
that the space is anisotropic but still homogeneous, so that there is no explicit dependence on the points x,x', but 
only on their difference. Namely, 



From the properties i) and ii) of Cap(r) one immediately obtains the following relations for the J-^s: 

Tt{r, z) ^ Ti{r,~z) 1,2,5 
Tsir, z) = -T^ir, -z) 
T'i{r,z) = Ti{r,z). 



(4.2) 



(4.3) 
(4.4) 
(4.5) 



Substituting the expression (4.2) into (3.14) and using the chain rules, we obtain, after lengthy but straightforward 



algebra, the following four equations (corresponding to the projections over raVp/r , 5ap, B^rp/r and B° B^): 

[air'^d^ + birdr + Ci(l - z^)9f + dizd^ + ei]J^i + 



[firdr + gizdz + ji] T2 + {kizrdr + hz^d^ 
[oi+piz'']n = {qi+riz')B°^ 



miz + nidz]J^3 



(4.6) 



02^1 + [b2r^d^ + c^rdr + ^2(1 - 7?)dl + 622^^ + h\ ^2 
92Z Tz + [fc2 + l2-^\ H = {m2 + n2Z^) B° ^ 



(4.7) 



[csr^d^ + d3rdr + 63(1 
[jaz rdr + ih + hz'')d, + m^z] T5 = n3B°^z 



z'')dl + hzd,+g3\T3 



(4.8) 



[biT^dl + CAvdr + ^4(1 — z'^)d1 + e^zdz + fi] J^5 — g^B 



(4.9) 

where the coefficients a^, bj , ■ ■ • are functions of ^ and d and are reported in Appendix Without loss of generality. 



we have fixed 13 = 1 in (2.4), and we have neglected all terms involving the magnetic diffusivity kq, our attention 
being indeed focused in the inertial range of scales, i.e., rj <^ r L where 77 = Kq^^ cx is the dissipative scale for 

the problem. 

With the substitution of the expression (4.2), the solenoidal condition ( [3.1?: ) splits into the following couple of 
equations: 



[ rdr + {d-l)]Ti + [rdr - zd,] T2 
\ z rdr + dz — z^dz — z\Tz — Q 



(4.10) 



dzT2 + \rdr + d\H + [z rdr + (1 - z')^.] .F5 = , 



(4.11) 



associated to the projections over rp/r and B °g, re spectively. 

From the relation (4.?;), Eqs. (4.10) and (4.11) it then follows that only two functions of the !F' s in ( ^^ ) are 
independent. A possible way to isolate contributions of the anisotropic components from the isotropic scaling is to 
use the decomposition of T' s on the Legendre polynomial basis. This is the subject of the next subsection. 
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A. Decomposition in Legendre polynomials 



In terms of the Legendre polynomials, functions J-i(r, z) can be decomposed in the form: 

oo 

T,{r, z)^Yl * = 1: 2, 5 (j even) 



^3(r,z)=^/f (r)F,(z) (j odd), 



(4.12) 
(4.13) 



where the separation of ev en a nd odd orders in ( 4.12 ) and ( 4.15 ) arises as a consequence of the symmetries expressed 
by the relations ( |4.3| ) and (4.4), respectively. 

Simple considerations related to the 'uniaxial' character of the forcing term with B° in the basic equation (^.l|), 
and the linearity of the latter in B° and B suggest that the index j in the above decompositions should be restricted 
to j < 2. The rigorous assessment of this point will be given in Subsec. IVB. On the other hand, contributions 



associated to j > 2 can be easily "activated" eithe r wh en a fully anisotropic forcing (i.e., projecting onto all Legendre 
polynomials) is added on the right hand side of ( |2.l|) , or in the framework of finite-size systems led by anisotropic 
boundary conditions. Moreover, as we shall see in Sec. ^ scaling exponents associated to j > 2 contribute to the 
inertial-range scaling of higher-order correlation functions involving the product BaBp at a single spacetime point. 
The latter property holds also without the invocation of a fully anisotropi c forc ing on the left hand side of Eq. (^.l|). 
From all these considerations we shall exploit the general decompositions ( 4.12| ) and (4.12) involving all j's. 



In order to obtain equations for /, 



we have to insert Eqs. (4.12) and (4.13) into Eqs. (|4.6|)-(4.11). Furthermore, 



quantities like z^d^T (p = 0, 1, 2 and q = 0, 1) and (1 — z^^d^T have to be expressed in terms of Legendre polynomials. 
This can be done exploiting well-known relations involving the Legendre polynomials (see, e.g., Ref. [Q): the resulting 
expressions for z^djJ- and (1 — z^)d^J-' are reported in Appendix M. For the sake of brevity, we report hereafter only 



the projection of equation (4.(;), the structure of the others being indeed similar (the full set of equations is however 
reported in Appendix O 



- Cl 



/(l-j)/«+2(2j + l)^/,«^. 



9=1 



./j^^ + (2. + l)Ei^^l. 



9=1 



ei/. 



(1) 



'(2) 



9i 



.4^^ + (2. + DEC 

9=1 



+ Jl/. 



:2) 



kir 



2j-l 



/(3) I J + 1 /(3) 

Jj-i ^ On I Q-'i+l 



Jjj - 1) .(3) (j + l)(j+2) ^(3) 
J,-l 



2j-l 



2r 



B 



J 

o 2 



(2j-l)(2j-3) 



A5) , 

Jj-2 + 



/j^\ + (2j + i)E4Vi 

9=0 



mi 



2j + 3 

i 



ii(2j-f l)E/2%+i 

9=0 



2j-l 



(2j + 3)(2j-|-l) (2j-l)(2j + l) 



.(3) I i + 1 .(3) 

/j-i + 2j-K3-^^'+\ 



.(5) (.? + 2)(j + l) (5)- 

^= ^ (2j + 3)(2j-f 5)^^+^. 



9i 



^1 1 3^J,2 



-^3.0 



(4.14) 



From the above equation we can see that terms like z^d^T are responsible for the coupling between an arbitrary 



anisotropic contribution of order j and all larger orders. The full set of equations (C1)-(C6) is thus not closed and there 
are no chances to solve them analytically. Simple physical argumentations actually permit to overcome the closure 
problem. Indeed, in the presence of a cascade-like mechanism of energy transfer towards small scales, anisotropy 
present at the integral scale should rapidly decay during the multiple-step transfer, and an almost isotropic inertial 
range scaling should be restored. Mathematically, this means that /j''''s should be rapidly decreasing functions of the 
degree of anisotropy j, i.e.. 



/« « !f « 



(4.15) 



and similarly for their derivatives. We shall control the validity of this physical assumption in a self-consistent way, 
at the end of our calculation.]^ 



^The physical assumption (4.15) is unnecessary when the decomposition in the irreducible representations of the SO{d) 
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The hierarchy (4.15) is exploited here by retaining, for each i appearing in the functions fp'^^s, the lowest value of 
the index p. When doing this, the simplifications on Eqs. (C1)-(|C6|) are enormous and the resulting set of equations 
reads: 



jU - 1) 



(1) 



2j-l 



+ M(l - j) + dij + ei] /j'^ + hrf}'^ + + ji] + kir^r^f'/^l 



2j-l 



A3) 



Pi 



(2j-l)(2j-3) 



A5) 
J3--2. 



B 



qi + ri 



2j- 



(4.16) 



(2j-l)(2j-3) 



.92 



(3) 



a3(2j + l)/«i 



2j 



a4(2j + l)/j^\ 



2j 



+ 



f(5) 



o 2 



m2 



'2 1 



.2f"(3) 



(3) 



p(3) 



J(J - 1) 
.2f"(5) 



TO3 



2j-l 



^(5) 



J,2 



(5) 



.94 B° '^,,0 



(4.17) 

(4.18) 
(4.19) 



^'(1) I 1^-f(l) , ^'(2) ■/(2) , .? --'(3) J f(3) ^ 



(2j + l)/j^\+r//^^+d/j^^ + 



'(3) 



f(3) 



2j 



2j 

j(j - 1) ^(5) 

2j 



1 ^-1 



0. 



(4.20) 
(4.21) 



Some remarks are worth. Focusing on the isotropic contribution, j = 0, we notice that the first two equations 
involve solely the functions /q^' and /q'^-' (and their derivatives). With the solenoidal condition ( 4.20 ), it is easy to 
check that they coincide with the equation reported in Ref. p7[ for the isotropic problem. Moreover, for j > 2, 

Eq. (4.16) suggests to take the function fj = {fj^\fj^\fj^i,fj^2) unknown field. It is immediately verified 

that fj appears al so in the other equations when the index j {j > 2) is renamed {j — 1) in Eqs. (4.18) and (4.21) 
and {j — 2) in Eq. (4.19). When doing this, a linear partial differential equation (PDE) system of the type Cjfj = Qj 
(hereafter, repeated indices are not summed) is obtained, Qj involving all terms related to the mean field B° , and Cj 
is restricted, for instance, to the first four equations. 

The analytical treatment of the resulting equation system still remains a very hard task for general values of the 
space separation r. The situation changes when one focuses on the inertial range of scales (i.e., for rj r L). In 
the latter case scaling behaviors are expected and we shall have: 



(r) oc r^j 



with 



(4.22) 



where the hierarchy on the exponents (j*'' immediately follows from (4.15). 

The structure of the above equations fixes the relation between the scaling exponents relative to different /'s. 

Indeed, when searching for power law solutions ff\r) oc r^j , in order to obtain balanced equations the 'oblique' 
relations must hold: 



^ c 



(1) 



(2) 



^(3) 



A5) 
^3-2- 



(4.23) 



We are now ready to show that nontrivial scaling behaviors for /, t ake place due to zero modes, i.e., solutions of the 
homogeneous problem Cjfj ~ 0. To that purpose, we exploit (4.23) and define coefficients yj through the relation 
fj = yjT^^ ■ Inserting the latter expression in the PDE system, a 4 x 4 algebraic linear system for yj is obtained. 



symmetry group is exploited (see Ref. |35|] for the case d = 3). This leads exactly to the same results obtained earlier in 
Ref. ||l^, where the hierarchy (4.15) was assumed. Notice that the additional exponents (subsets II and III) reported in ^] 
are related to the pseudotensorial structures and that in our model they do not contribute to the inertial-range behavior of the 
pair correlator. 
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The emergence of zero modes is thus reduced to impose the existence of nontrivial solutions of a 4 x 4 homogeneous 
Unear system, that means here the resohition of an algebraic equation of 8-th degree arising from the condition that 
the determinant of the system coefficients is zero. The calculation, lengthy but straightforward, leads to four sets of 
zero modes (actually eight s ets, b ut it turns out that the associated coefficients, y^, of four of them do not satisfy the 
solenoidal condition ( [4.2C| )-( f.21 )) the expressions and the admissibility of which are given and discussed in the next 
section. 



B. Zero-mode solutions and their admissibility 



Let us start with the case j = 0, corresponding to the isotropic contribution. As already observed, Eqs. (4.16) and 
( 4.17 ) are decoupled from the others, and the problem can be solved directly for /g^"* and /g^"* which must satisfy 



also the solenoidal condition (4.20). The imposition of the existence of nontri vial solutions (for j ~ we have a 
homogeneous 2x2 algebraic linear system for yj ) and the solenoidal constraint ( 4.20|) lead to the following solutions: 



with their ^ — > and d — > oo limits: 



2d - 2 



Co+--e + o(e^) 
= -e-^ + o(i/d^) 



d- 3 



d- 1 

d d 



o{e) 

^- + 0(l/d^). 



For j > 2 zero-mode exponents are: 



1 



d^ -d- 



2(d 

6d^ +Ai^d + 8 + 10di + 20dj 



-2d^£_-2d^i^ 



20d- 8C- 
j + Ad^ f + 2^2 „ 4^j2 _^ -^7^2 „ 8dj2 + 

ii + Ad^ j + Ad^ ii + Adf i + Af ~ IQd^j 

, -|l/2~ 

nd£,j + d'^± 2Vk {d - I) (2 - £,) 



(4.24) 



(4.25) 



(4.26) 



(4.27) 



where 
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d^-d + 



- 2 d^ ^ - 2 d^ ^2 



2(d 

6d^+Ai^d + 8 + 10di + 20dj-20d~8i- 
8j + Ad'^f + 2^^ -A^f + nd^ ~Sdf + 
S^j + Ad^j + Ad'^jC + AdfC + '^f - 16d2j 

, -|l/2~ 

ndU + d'^± 2VZ(d- 1) (2-e) 



(4.28) 



K = {d-l){d^ + Ad^ j -5d^ + 2d^i + i^ d + AdU 



6di + 8d-l2dj + Adf 
4-4/ -f 4^-2) 



i^+A^ + Sj-Sa 



with their ^ — > and d — > oo limits: 
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J s 

3-i 



{d-2 + 2j){d + 2j){d- 1) 



„ ^ -4rf2 - i7dj + I6d + 28j - 16 - Uf + 2d'^j + 5df + 2f 
-J -2 + ^ ^-^^ — ^ ^ ^ ^ , ^ . — ^^^7^ — ^ — + 0{i^) 



2^3 



1 



(d-2 + 2j)(d + 2j-4)(d-l) 
+ 0(1/^2) 

-5^2 - 7dj + 6d + Aj - 2p + d^ + 2d'^j - df - 2f 



id-2 + 2j){d + 2j){d~l) 
'-^^l±^ + Oil/d^) 



o{e 



p, = 



^ , , rj-l)(2j2-4j + 5dj-4d+2d2) 



1 



+ 2 -J 



(d-2 + 2j)(rf + 2j-4)(d-l) 
'^^'~'K0il/d^). 



(4.30) 



(4.31) 



(4.32) 



d d 

Let us discuss the admissibility of these solutions. With the term admissible we mean a solution fj{r) satisfying 
the appropriate boundary conditions, at both small (UV limit) and large scales (IR limit). Specifically, the following 
asymptotic behaviors have to be satisfied: 



fj{r) regular for r ry — > 







for 



r » L. 



(4.33) 
(4.34) 



Concerning the limit (4.33), we have to consid er so lutions corresponding to the diffusive range and to match them 
with our inertial range power laws. From Eq. (3.14) we can easily see that the equations holding in the diffusive 
range are obtained by setting to zero the parameter ^. The consequence is that our inertial range zero-mode solutions 
become solutions in the diffusive range for ^ = 0. The problem related to the UV boundary condition is thus reduced 
to search for regular solutions for fj{r) in the ^ ^ limit. This is easy to do and the result is that solely exponents 
(^Q and for j > 2 permit to satisfy the condition of regularity for fj, the other exponents being indeed < for 

^ = 0. Notice that, the zero-mode exponent (o coinci des w ith the isotropic solution obtained in Ref. p7| ]. 

Let us pass to discuss the IR boundary conditions ( 4.34 ). In this case, as pointed out in Ref. a crucial role is 
played by the external forcing. Indeed, in the presence of forcing, zero modes and the decaying forced solution may be 
matched at the integral scale L, thus satisfying the IR boundary conditions. The result of this argument (which can 
be rigorously illustrated solely for j = where the general solut ion fo r /, (r) is available) is that zero-mode exponents 



are not admissible for j > 4. Indeed, as we can see from Eq. (4.16)-(4.19), the forcing term related to B" projects 
solely on the shells j < 2. 

To summarize, we have one admissible zero mode for j — (Co) ^-nd two admissible zero modes for j — 2 (C^). 
Our attention being focused on the inertial range of scales (i.e., r/L << 1), our choice for j — 2 is for ■ We have 
indeed to take the exponent giving the leading inertial range contribution. 

Finally, we can thus define the final solution Q of our problem as: 



Co ^ C = 



+ d - 2^ + Vl2d2^ - Sd^ + S^d - M^^"^ - Ad^( + d"^ + d^ - 2d^ 



2d -2 



= -^-^ + 0{l/d') 



(4.35) 



C2 ^ Co 



2{d- 
2^ + 2d^ + 
2^ 

(d- l)(d + 2) 
§ + 0(l/d^) 



d + 2^ 

--2\/x(2-C)(d-l) 



8 - 12d - 8C + 2di + Aifd - 2d^£,^ - 2d^^ + S^C + 
1/2' 



(4.36) 
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where 



K = {d- l){d^ + 2(f£, + + + 2rfe - 4 + 4^ - i^). 



We stress that, for j > 2, exponents Cj = C/ become admissible under the conditions aheady discussed in Subsec. IV A. 

The last remark concerns the self-consistency of our solution, that is, the validity of the hierarchy in ( 4.22 ). The 
validity of the latter can be easily verified from Fig. |l| where the behavior of C,j {j = 0, 2, 4 and 6) are shown for d = 3 



as a function of ^. Similar behaviors actually hold for all d's and j's. As we shall see in Sec. |V|, a hierarchical order 
for the scaling exponents is also present for higher-order correlation functions. 



V. UV RENORMALIZATION OF THE MODEL. RG FUNCTIONS AND RG EQUATIONS 

The RG approach to the statistical models of the turbulence is exposed in Refs. p^ , ^ in detail (see also Ref. p^ ] 
for the scalar Kraichnan model); below wc confine ourselves to the only information we need. 

The analysis of UV divergences is based on the analysis of canonical dimensions (see, e.g., Ref. Isi)). Dynamical 



models of the type (3.1), in contrast to static models, have two scales, i.e., the canonical ("engineering") dimension 
of some quantity F (a field or a parameter in the action functional) is described by two numbers, the momentum 
dimension and the frequency dimension df^. They are determined so that [F] ^ [£]~'^p [T] where i is the length 
scale and T is the time scale. The dimensions are found from the obvious normalization conditions 

4- -4-1, d- = < = o, dt^d^^o, d^ = -dt = i, 

and from the requirement that each term of the action functional be dimensionless (with respect to the momentum 
and frequency dimensions separately). Then, based on dp and dp, one can introduce the total canonical dimension 
dp = dp + 2dp (in the free theory, dt oc d^), which plays in the theory of renormalization of dynamical models the 
same role as the conventional (momentum) dimension does in static problems. 



In the action (3.1), there are fewer terms than fields and parameters, and the canonical dimensions are not deter- 



mined unambiguously. This is of course a manifestation of the fact that the "superfluous" parameter B° = \B°\ can 



be scaled out from the action (see Sec. III). After it has been eliminated (or, equivalently, zero canonical dimensions 
have been assigned to it), the definite canonical dimensions can be assigned to the other quantities. They are given in 
Table |, including the dimensions of renormalized parameters, which will appear later on. From Table | it follows that 
the model becomes logarithmic (the coupling constant go becomes dimensionless) at ^ = 0, and the UV divergences 
have the form of the poles in ^ in the Green functions. 

The total canonical dimension of an arbitrary 1-irreducible Green function F = ($ • • • <l')i-ir is given by the relation 

dr = d^, + 2d^ = d + 2- N^di, (5.1) 

where 7V$ = {Nb' , Nb, N^} are the numbers of corresponding fields <& = {B',B,v} entering into the function F, 
and the summation over all types of the fields is implied. The total dimension dr is the formal index of the UV 
divergence. This means that superficial UV divergences, whose removal requires counterterms, can be present only in 
those functions F for which dr is a non-negative integer. 

Analysis of divergences in the problem (|3.1|) should be based on the following auxiliary considerations: 



(i) All the 1-irreducible Green functions with Nb> < Nb vanish (see Sec. III). 

(ii) If for some reason a number of external momenta occur as an overall factor in all the diagrams of a given Green 
function, the real index of divergence d'^ is smaller than dr by the corresponding number (the Green function requires 
counterterms onl y if dp is a non-negative integer). 

In the model ( |3.lD , the derivative d at the vertex can be moved onto the field B' using the integration by parts, 
which decreases the real index of divergence: dp = dr — Nb' ■ The field B' enters into the counterterms only in the 
form of a derivative, daB'^. 



(iii) A great deal of diagrams in the model ( |3.l[) contain effectively closed circuits of retarded propagators {BB')o 
and therefore vanish. For example, all the nontrivial diagrams of the 1-irreducible function {BaB'^v.y)i-ir vanish. 

From the dimensions in Tabic | we find dr = d -|- 2 - 7V„ - dTVs' and d^ = (d -I- 2) - A'^y - (d -I- 1)Nb'. From 
these expressions it follows that for any d, superficial divergences can only exist in the 1-irreducible functions with 
Nb' = 1,N,=Nb = (dr = 2, d'^ = 1), Nb' =Nb = 1,N,=0 (dr = 2, dj, = 1), Nb' =N, = 1,Nb=0 (dr = 1, 
dp = 0), and Nb' = Nb = N,, — 1, (dr = 1, dp — 0) [we recall that Nb < Nb' (see (i) above)]. However, the first of 
these counterterms has necessarily the form of a total derivative, B'^d^Ba, vanishes after the integration over x and 
therefore gives no contribution to the renormalized action. Furthermore, for the last two of these functions, all the 
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nontrivial diagrams vanish [see (iii) above]. As in the case of the passive scalar field |2^, we are left with the only 
superficially divergent function the corresponding counterterm necessarily contains the derivative d and 

therefore reduces to B'^d'^Ba (another structure, B'^dadjsB/^, vanishes by virtue of the solenoidality of B). 

Introduction of this counterterm is reproduced by the multiplicative renormalization of the parameters ga, kq in the 
action functional ([3.1|) with the only independent renormalization constant Z^: 



Kq 



kZ, 



K, .go = ffM 



z„ = z: 



(5.2) 



Here fi is the reference mass in the minimal subtraction scheme (MS), which we always use in what follows, g and 
K are renormalized analogs of the bare parameters go and kq, and Z — Z{g,£^,d) are the renormalization constants. 
Their relation in Eq. (^.2| ) results from the absence of renormalization of the contribution with K^^ in Eq. (3.1), so 
that Do = goKQ = giJ^K. No renormalization of the fields and the "mass" m is required, i.e., Z$ = 1 for all 
mo = m. The renormalized action functional has the form 



$ and 



5i?($) = B' -dtB - [v ■ d)B + {B ■ d)v + {B° ■ d)v + Z^nd^B - 



(5.3) 



where the function K from Eq. ( 2.21 ) is expressed in renormalized parameters using Eqs. (5.2): _Do = .go'^o ~ SM^'*- 
The relation S'($, ep) = 5'ij(<i>, e, /i) (where eo = {go, kq, m} is the complete set of bare parameters, and e = {g, k, m} 
is the set of renormalized parameters) implies W{eo) = WR{e, /i) for the bare correlation functions W = ($•••$) and 
their renormalized analogs Wr. We use I?^ to denote the differential operation for fixed eg and operate on both 
sides of this equation with it. This gives the basic RG differential equation: 



(5.4) 



where Vrg is the operation I?^ expressed in the renormalized variables: 

VRG = V^+f3{g)dg--f^{g)V^. (5.5) 

In Eq. (|5.5| ), we have written = xdx for any variable x, and the RG functions (the /3 function and the anomalous 
dimension 7) are defined as 



"fpig) = 'D^^'o^Zp for any Zp, 



(5.6a) 



Pig) = V^g^g[-^ + j^ig)]. 



(5.6b) 



The relation between /3 and 7 in Eq. ( 5.6b| ) results from the definitions and the last relation in Eq. ( ^.2[ ). 

Now let us turn to the explicit calculation of the constant Z^ in the one-loop approximation in the MS scheme. It is 
determined by the requirement that the 1-irreducible function expressed in renormalized variables be UV 

finite (i.e., be finite for ^ — > 0). This requirement determines up to an UV finite contribution; the latter is fixed 
by the choice of a renormalization scheme. In the MS scheme all renormalizatio n con stants have the form "1 + only 
poles in The function { B'B) i^i^ in our model is known exactly [see Eqs. ( 3.4a ) and ( 3.11 )]. Let us substitute 
Eqs. (5.2) into Eqs. (3.4a), (3.11) and choose Z„ to cancel the pole in ^ in the integral J{m). This gives: 



Z^ = l-gC, 



2d^ 



(5.7) 



5^ 



with the coefficient Cd from Eq. ( |3.11b| ). Note that the result ( ^.7[ ) is exact, i.e., it has no corrections of order 
g^, and so on; this is a co nseq uence of the fact that the one-loop approximation (3.11) for the response function is 
exact. Note also that Eq. (5.7) coincides literally with the exact expression for Z„ in the case of a passive scalar (see 
Ref. 

For the anomalous dimension 7k(<?) = T^^i^^Z^ — j3{g)dg\nZ^ from the relations ( ^.6bD and (5.7) one obtains: 



7«(.9) 



gCd 



jd-l) 
2d 



(5.8) 



From Eq. (5.6b) it then follows that the RG equations of the model have an IR stable fixed point [(3{g*) — 0, 
P'{9*) > 0] with the coordinate 
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Let F(r) be some equal-time two-point quantity, for example, the pair correlation function of the primary fields 
<f> = {-B', B, v} or some composite operators. We assume that F{r) is multiplicatively renormalizable, i.e., F = ZpF^ 
with certain renormalization constant Zp. The existence of nontrivial IR stable fixed point implies that in the IR 
asymptotic region Ar >> 1 and any fixed mr the function F(r) takes on the form 

F{r) ~ nf^ A'^'-iAr)-^^ ximr), (5.10) 

where dl^ and dp are the frequency and total canonical dimensions of F, respectively, and x is some function whose 
explicit form is not determined by the RG equation itself. The critical dimension Aj? is given by the expression 

A[F] = Ap = d'p + A^d^p+-f*p, (5.11) 



where jp is the value of the anomalous dimension (5.6a) at the fixed point and A„ = 2 — 7* = 2 — ^ is the critical 
dimension of frequency [note that the value of 7^(5) at the fixed point is also found exactly from the last relation in 
Eq. I^y. 7: = 7«(5,)=^]. 



The critical dimensions of the basic fields in our model are found exactly [we recall that they are not renormalized 
and therefore 7$ = for all $]. From the dimensions in Table | we then find A„ = 1 — ^, Ag = 0, Ag/ = d. 

VI. CRITICAL DIMENSIONS OF COMPOSITE OPERATORS 

Any local (unless stated to be otherwise) monomial or polynomial constructed of primary fields and their derivatives 
at a single spacetime point x = {t,x} is termed a composite operator. Examples are B^, B'^d^Bp, B'^v ■ dBa, and 
so on. 

Since the arguments of the fields coincide, correlation functions with these operators contain additional UV di- 
vergences, which are removed by additional renormalization procedure (see, e.g., Ref. [0). For the renormalized 
correlation functions standard RG equations are obtained, which describe IR scaling with definite critical dimensions 
Ap = A[F] of certain "basis" operators F. Owing to the renormalization, A[F] does not coincide in general with the 
naive sum of critical dimensions of the fields and derivatives entering into F. 

In general, composite operators are mixed in renormalization, i.e., an UV finite renormalized operator F^ has the 
form F^ — F+ counterterms, where the contribution of the counterterms is a linear combination of F itself and, 
possibly, other unrenormalized operators which "admix" to F in renormalization. 

Let F = {Fa} be a closed set, all of whose monomials mix only with each other in renormalization. The renormal- 
ization matrix Zp = {Zat} and the matrix of anomalous dimensions "fp = {jab} for this set are given by 

Fa = J2 ^-bF^^ TF = Zp^V^Zp, (6.1) 



and the corresponding matrix of critical dimensions Ap = {Aab} is given by Eq. ( |3.11 ), in which dp, df,, and dp are 



understood as the diagonal matrices of canonical dimensions of the operators in question (with the diagonal elements 
equa l to sums of corresp ond ing dimensions of all fields and derivatives constituting F) and j*p = jp{g*) is the matrix 



(3T) at the fixed point (5.9) 



Critical dimensions of the set F = {Fa} are given by the eigenvalues of the matrix Ap. The "basis" operators that 
possess definite critical dimensions have the form 

F^^^ ^Y.UabFi' , (6.2) 

b 

where the matrix Up — {Uab} is such that A^ = UpApUp^ is diagonal. 

In general, counterterms to a given operator F are determined by all possible 1-irreducible Green functions with 
one operator F and arbitrary number of primary fields, F — {F{x)^{xi) ■ ■ ■ ^{xn))i~ir- The total canonical dimension 
(formal index of divergence) for such functions is given by 

dr = dp — Af$c?$, (6.3) 

with the summation over all types of fields $ = {B',B,v} entering into the function. For superficially divergent 
diagrams, the real index of divergence, dp = dp — Nb', is a non-negative integer, cf. Sec. 0. 
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In what follows, an important role will be played by the tensor composite operators built solely of the field B 
without derivatives: 



-P- 



(6.4) 



Here p is the rank of the tensor and n = 21 + p is the total number of fields B entering into the operator. From Table 



I and Eq. ( |6.3| ) for the operators (3.4) we obtain dp = and dp = —N^ — {d + 1)Nb'- Therefore, the divergences 



can exist only in the functions wit h N ^, = Nb' = 0, for which dr = d'y = This means that the operators F^^p'^ mix 
only with each other, i.e., the set (6.4) is closed with respect to the renormalization. 
The simple analysis of the diagrams shows that the 1-irreducible function 



(F("f)(x)B(xi)---B(.T„0)] 



(6.5) 



contains the factor (i?°)" " and therefore vanishes for n' > n, 
follows that the operator F^^" p '> can admix to F'"^^ only if n' 
renormalization matrix 



Cf. 

< n 

(np) 



the discussion in the end of Sec. III. It then 



This means that the corresponding infinite 



(6.6) 



n' p' 



is in fact block-triangular, i.e., Znp^n'p' — for n' > n, and so are the matrices 7^?, A^? and Up- It is then obvious that 
the critical dimensions associated with the operators F^"'') are completely determined by the eigenvalues of the finite 
subblocks with n' = n. In the following, we shall not be interested in the precise form of the basis operators (|6.2|), 
we rather s hall be interested in the anomalous dimensions themselves. Theref ore , we can neglect all the elements of 
the matrix (^^) other than Znp,np'- The latter are found from the functions which are independent of B° and 
therefore can be calculated directly in the isotropic theory with B° — 0. It is then clear that the block Znp^np' can be 
diagonahzed by the changing to irreducible operators: scalars {p ~ 0), vectors {p = 1) and traceless tensors {p > 2), 
but for our purposes it is sufficient to note that the elements of the block Z^p^np' vanish for p < p' , i.e., this block 
is triangular along with the corresponding blocks of the matrices 7^^, Aj? and Up- Indeed, the irreducible tensor of 
the rank p consists of the monomials with p' < p only, for example, F^^^-* = BaBfs — SapB'^/d, and therefore only 
these monomials can admix to the m onom ial of the rank p in renormalization. The final conclusion is that the critical 
dimensions, associated with the set ( |6.4| ), coincide with the diagonal elements An,p = Anp^ np of the matrix (5.11), 



they are completely determined by the diagonal elements Znp = Znp^np of the matrix (3.6), and that they can be 
calculated directly in the isotropic theory with B° = 0. 

Now let us turn to the one-loop calculation of the diagonal element Z„p of the matrix Zp in the MS scheme. Let 
r(a;; B) be the generating functional of the 1-irreducible Green functions with one composite operator F'^"^'^ and 
any number of fields B. Here x = {t, x} is the argument of the operator and B{x) is the functional argument, the 
"classical analog" of the random field B{x). We are interested in the n-th term of the expansion of r(a;; B) in B{x), 
which we denote T^'^'>{x; B); it has the form 



Toi ■■■Qp (2^1 -B) 



- y 



dxi 



dXnBf3^{xi) ■ ■ -Bp^iXn) {F^al^lap{x)B I3^{xi) ■ • •B^^(x„))i_ 



(6.7) 



In the one-loop approximation the functional (6/7) is represented diagrammatically as follows: 



F 



(np) 



(6.8) 



Here the solid lines denote the hare propagators {BB)q from Eq. (3.3), the ends with a slash corres pon d to the field 
B' , and the ends without a s lash correspond to B; th e da shed line denotes the velocity propagator (2.2); the vertices 
correspond to the factor (3^). The first term in Eq. ( |6.8| ) is the "tree" approximation, and the black circle with two 
attached lines in the diagram denotes the variational derivative 



14i---ap/3i/32(2;; Xi,X2) 



^ -^Oii ■ ••Q.p (*^) 

5Bp^[xi)5Bp^ (X2)' 



(6.9) 



It is convenient to represent it in the form 
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ipftftC^;; xi,x2) = S{x - xi)d{x - X2) 



dbfj, db 



132 



bc.,---bc.^ {b'^y 



(6.10) 



where ba is a const ant vector, which after the differentiation is substituted with the field Ba{x). 

The vertex ( 6.1C ) contains (n — 2) factors of B. Two remaining "tails" B are attached to the lower vertices of the 



diagram in Eq. (6.8). We know that the UV divergent part of the diagram is proportional to n factors B without 
derivati ves^ so that we can omit the first term of the vertex B' [—{v ■ d)B + {B ■ d)v\^ or, equivalently, the first term 
in Eq. ( |3.2| ). Further mor e, we can set all the external momenta in the integrand equal to zero, and the UV divergent 
part of the diagram (|6.8[) takes on the form 



bp-^bp^ 



92 



dbn, db, 



where we have denoted 



T, 



1310213304 



•Pl ^"P2 



Do 



T, 



/3i/32/33/34 i 



duj 
2^ 



qp3<ii3i PpiP2iq) 



(6.11) 



(6.12) 



We recall that the integration over q should be cut off from below at q = m (see Sec. In Eq. ( 6.12 ), we have to 
change to the renormalized variables using Eqs. (5.2); in our approximation this r educes to the substitution go ~^ 91^-^ 
and kq — > K. Then we perform the integration over lj and use the relations ( |3.1C| ) and 

^gftgfegfeg/?4 _ ^l3ll32^Pif3i + ^PlPi^P2[33 + hip3^l32l3i 



dqfiq)- 



d{d + 2) 



dqfiq). 



This gives 



T, 



/3i/32/93/34 



g^^ J(m) 
2d{d + 2) 



{d + l)6p^pjp^p^ - {5pipJp2P3 + Sp^pjp^p^) 



(6.13) 



with the integral J{ m) de fined in E q. (|3.11b|) . 

Substituting Eq. (|6T|) into Eq. ( |6.1lD gives the desired expression for the divergent part of the diagram (p^). It 
is sufficient to take into account only the terms proportional to the monomial B^^ (x) ■ ■ ■ Ba {x) [Ba{x)Ba{x)Y and 
neglect all the other terms, namely, those containing the factors of 5aia2 stc. The latter determine non-diagonal 
elements of the matrix Zp, which we are not interested in here. Finally we obtain 



(np) 



gfi^ J{m) 



'np 



M{d + 2) 



(6.14) 



where we have written 



2n{n-l)- {d+l){n-p){d + n+p-2) 
2p (p - 1) - (d - 1) (n - p) (d + n + p). 



(6.15) 



The dots in Eq. (6.14) stand for the 0{g^) terms and the structures different from F^"p\ ~ denotes the equality up 

to UV finite parts; we also recall that n = p + 21. 

The constant Znp is found from the requirement that the renormalized analo g F^ = Z^^Tn of the function (6.14) 
be UV finite (mind the minus sign in the exponent); along with the expression ( 3.11b ) for the integral J(m) and the 
MS scheme this gives 



Qnp 

4d{d + 2)^ 



0(9' 



with Cd from Eq. (3.11b). For the anomalous dimension j„p — T^ii ^T^Znp it then follows 

.gCd Qnp 



Inpig) 



Ad{d + 2) 



0{9'). 



(6.16) 



(6.17) 



From Table | and Eqs. (5^) and ( 5.11 ) for the corresponding critical d iniension p — ^np(y*) finRlly obtciin. 



16 



(6.18) 



with the polynomial Qnp from Eq. ( |6.15 ). 

The straightforward analysis of the expression ( |6.1£ ) shows that for fixed n and any d, the dimension A„^p 
monotonically with p and reaches its minimum for the minimal possible value, i.e., p = if n is even and p 
is odd: 



if p > p . 

Furthermore, this minimal value is negative and it decreases monotonically as n increases: 

> A2fc,0 > A2i:+14 > A2fc+2,0- 



Finally, we note that for any fixed p, the dimension ( 3.1q) decreases monotonically as n increases: 



if n < n' 



decreases 
1 if n 



(6.19a) 
(6.19b) 
(6.19c) 



The inequalities (6.19) show that the critical dimensions of the tensor operators ( |6.4| ) exhibit a kind of hie rarchy; in 
particular, the less is the rank, the more negative is the dimension and, as will be explained in Sec. VD, the more 
important is its contribution to the inertial-range behavior. 



In the model of passive scalar advection by the rapid-change velocity field (2^) in the presence of an imposed 
linear gradient, similar inequalities are satisfied by the critical dimensions of tensor operators of the type (6.4), but 
constr ucte d of gradients of the scalar field (see Ref. |l4j). In the order O(^) their critical dimensions coincide exactly 
with ( ^.1^ ), which is, however, an artifact of the one-loop approximation (see Sec. |VI 

As already said above, the operators that possess definite critical dimensions ( 
the basis operators related to the latter by the relations (6.1) and 



M)- 



T8D are not (3.4) themselves, but 
In the isotropic case {B° = 0), the basis 



operator with the dimension A„^p is a p-th rank traceless tensor constructed of all the monomials i^^" ^ ^ with n' = n 
and p' < p. When the background field B° is "turned on," the admixture of the monomials with n' < n and p' > p 
becomes possible. The "missing" fields B in the monomials with n' < n are substituted with the constant fields B° 



(the total number of the fields B and B° has to be equal n, owing to the linearity of the basic equation (2.1) in B 
and B°), while the "superfiuous" indices of the monomials with p' > p are contracted with the indices of B° , so that 
the basis operator remains a p-th rank traceless ten sor. And vice versa, the unrenormalized monomial from 
( 3^ is a linear combination of the basis operators (|6.2[ ) with respective dimensions A„' ^p/ . The hierarchy relations 
( 3.19[ ) then show that the minimal dimension entering into is A„^p^, where p„ is the minimal possible value of 

p for a given n, i.e., p„ = if n is even and p„ = 1 if n is odd. 



VII. OPERATOR PRODUCT EXPANSION AND THE ANOMALOUS SCALING FOR THE 

CORRELATION FUNCTIONS 



The representation ( 5.10D for any scaling function x('Tir) describes the behavior of the Green function for Ar >> 1 
and any fixed value of mr. The inertial range corresponds to the additional condition that mr << 1. The form of the 
function x('7ir) is not determined by the RG equations themselves; in the theory of critical phenomena, its behavior 
for mr ^ is studied using the well-known Wilson operator produc t expa nsion (OPE); see, e.g., Ref. This 
technique is also applicable to the theory of turbulence; see, e.g., Ref. 



32 3 



According to the OPE, the equal-time product Fi(x)F2{x') of two renormalized composite operators at cc = (a; 
£c')/2 — const and r = a; — a;' has the representation 



Fi(x)i^2(a;') 



(7.1) 



where the functions Ca are the Wilson coefficients regular in and Fa are, in general, all possible renormalized 
local composite operators allowed by symmetry; more precisely, the operators entering into the OPE are those which 
appear in the corresponding Taylor expansions, and also all possible operators that admix to them in renormalization. 
If these operators have additional vector indices, they are contracted with the corresponding indices of the coefficients 

Ca. 



Without loss of generality it can be assumed that the expansion in Eq. (7.1) is made in basis operators (5.2) with 
definite critical dimensions A^. The renormalized correlation function {Fi{x)F2{x')) is obtained by averaging Eq. 
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( [7.l[ ) with the weight expSji, the quantities {Fa) appear on the right hand side. Their asymptotic behavior for m ^ 
i s fo und from the corresponding RG equations and has the form (Fa) oc m'^". From the operator prod uct ex pansion 
( [7.l[ ) we therefore find the foUowing expression for the scahng function xl™*") ™ the representation (5.1C) for the 
correlation hmction {Fi{x)F2{x')): 



X{mr) 



(7.2) 



where the coefficients Aa = Aa{mr) are regular in (mr)'^; they depend on ^, d and, in our case, on the cosine 
z = cos9 = B ■ r/r. 

Consider for definiteness the equal-time pair correlation function of the operators ( |6.4| ); their vector indices will be 
omitted in order to simplify the notation. For the leading term in the asymptotic region Ar >> 1 from the general 
expression (5.10) we obtain 



(F("'')(xi)F("'p''(x2)) = (Ar) 



Xnp,n'p' {tTIt)^ 



(7.3) 



with the dimensions An,p from Eq. ( |6.1S ) and certain functions Ykp.k'p' (wr ). We recall that the monomial (6.4) is a 
linear combination of basis operators possessing definite critical dimensions ( |6.18 ) with different values of the indices; 
we also reca ll that p„ is the minimal possible value of p for a given n, i.e., p„ = if n is even and p„ = 1 if n is odd. 
In Eq. (^^), only the leading contribution is displayed which is determined by the minimal dimensions entering into 
the operators on the left hand side (see the discussion in th e en d of Sec. VI). 

The leading term of the Taylor expansion for the function ( |7.3| ) involves the operators F^'^') from ( |6^ ) with k — n+n' 
and I < p + p'] higher-order terms involve tensors of arbitrary rank, built of the field B and its derivatives. The 
decomposition in renormalized operators gives rise to all the tensors F'*^'^ with k < n + n' and all possible values of p; 
the tensors with I > p + p' appear owing to the renormalization of the hi ghe r-order terms with derivatives. Therefore, 
the desired asymptotic expression for the function Xnp,n'p'{''Tir) in Eq. ( [7.3| ) in the region mr << 1 has the form 



Xnp, 



jp'irnr) = ^ ^ Afcj (mr)^*^'' -|~ 



(7.4) 



fc=0 



where A^j are coefficients dependent only on ^, d and z = cosO, and the second summation runs over all values of j, 
allowed for a given k. 

Some remarks are now in order. 

The leading term of the inertial -rcingG bdis-vior {jiit <C<C 1) of the function x.np,n'p' 

(mr) is obviously given by the 

contribution with the minimal dimension A^j entering into Eq. ([7.4|). 



The dots in Eq. (7^) stand for the contributions of the order (rnr)^+'^(^^ and higher, which arise from the senior 
operators, for example, Bd^B and so on. 

The operators F^^^^ with k > n + n' (whose contributions would be m ore important) do not appear in Eq. ([7.4|), 
because they do not appear in the Taylor expansion of the function (7.3) and do not admix in renormalization to the 
terms of the Taylor expansion. In other words, the number of the fields B in the operator Fa entering into the right 



hand sides of the expa nsions (7.1) can never exceed the total number of the fields B in their left hand sides. 

The expansion ( [7.4| ) is consistent with the Legendre polynomial decomposition of the type (4.2) or, in general, with 
the decomposition in irreducible representations of the rotation group, employed in Refs. ]12| , |13| , ^ . This becomes 
especially clear if the left hand side of Eq. (7.1) involves only scalar quantities. Then all vector indices of the mean 
values {Fa) in the right hand side are contracted with the indices of the corresponding Wilson coefficients Ca(r). As is 
explained in Sec. VI, the basis operator that possesses definite critical dimension Akj is a j-ih. rank traceless tensor, 
so that its mean value is also a j-ih rank traceless tensor, built solely of the constant vector B° and Kronecker delta 
symbols. It is then clear that its contraction with Ca(r) gives rise to the j-th order Legendre polynomial Pj{z). 

Now let us turn to the comparison of the nonperturbative results for the pair correlation function Cafsir) = (BaBfj), 
obtained in Sec. |^ using the zero-mode t ech niques, with the predictions of the RG and OP E, given above. To this 
end, we put n = rv~ p = p' = 1 in Eqs. ( [7.3| ), ([7.4]). The isotropic shell (j = 0) in Eq. (7^) is then represented by 
the tr ivial opera tor F = 1 (fc = 0) with Ao,o = and the monomial B^ = B^Ba {k — 2) with A2,o = + 0{^^) (see 
Eqs. ( 6.15 ) and ( 3.18 )). The leading term of the small-mr behavior is given by the latter, so that we have to identify 
A2,o with Co = ((T^om Eq. (|!p). 

It was mentioned in Sec. VI that in the one- loop approximation, dimensions ( |6.18 ) coincide with the critical 
dimensions of tensor operators of the type (6.4), but constructed of the scalar gradients. The above identification 



shows that this coincidence is confined to the order O(^) even for the simplest dimension A2,o- For the scalar case. 
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one has A2.0 = — ^ exactly, in agreement with the well-known exact solution for the two-point structure function 
obtained in |20[| , while in our case A2,o is a nontrivial function of ^. 

At first sight, the first anisotropic correction is related to the term with k = j = 1 in Eq. (7.4), i. e., to the simplest 
operator B. However, the mean value {B{x)) vanishes and therefore gives no contribution to Eq. (7.4). Indeed, the 
analysis of the diagrams shows that ( B(x )) is obtained from the 1-irreduciblc function {B' {x))i-n-, which vanishes 
owing to the invariance of the model {"iA) with respect to the shift B' ^ B' + const. 

The leading anisotropic correction is therefore related to the term with k = j — 2, i.e., with the operator BaJBg. 
Its dimension A2,2 = 2^/ (d —l){d-\-2) + 0(^^) has to be identified with ^2 = C2' ^'^'i agreement with Eq. ( 4.36 ). 

We have thus established the agreement between the O(^) results obtained using the RG and OPE, with the first 
terms of the expansions in ^ of the exact nonperturbative results obtained within the zero-mode techniques. Note 
that for the isotropic exp onent , such agreement was mentioned earlier in Rcf. to the order O(C^). 

The exact expressions (4.35), (4.36) can therefore be viewed as nonperturbative predictions for the critical dimen- 
sions of the operators B^ = BaBa and BaBp, respectively. Similarly, the results (4.27) for the higher exponents 
can be linked to certain composite operators with two fields B and j derivatives for (j' and {j — 2) derivatives for 
We shall not dwell on this point here and only note that the exponents (2 ^^.d are indeed related to the 
second-rank and fourth-rank families of the irreducible operators built of two fields B and two derivatives, dBdB, 
with various arrangements of the vector indices. 

As is explained above in Sec. [V, the exponents (^^ for j > 4 do not appear in the inertial-range behavior of the 
pair correlation function. This is also easily understood within the OPE. The mean value of the j-th rank irreducible 
operator with n fields S is a traceless j-th rank tensor built of n vectors B° and Kronecker delta symbols. Th is 
follows from the linearity of the basic equation ( ^.1[ ) in B and B° (see also the discussion in Sec. Ill below Eqs. ( |3.3[ )). 
However, nonvanishing tensors of this type do not exist if the number of vector indices exceeds the number of fields 
[the structures like B°^B°^B°^B°^/{B°)'^ are forbidden because B° appears in the bare propagators (3.3) only in the 
numerators] . 

It was noted in Sec. ^ that these exponents will be activated when a f ully anisotropic forcing term (i.e., projecting 



onto all Legendre polynomials) is added to the right hand side of Eq. (2.1). Mor eove r, the above interpretation in 



terms of the OPE shows that they are relevant even for the original simple model (2.1). 

A lthough the contributions with j > n vanish in the mean value {Ba(x)Bp{x')), they are present in the expansion 
(7.1) without averaging and therefore the exponents (j*^ can reveal themselves in other correlation functions that 
involve the product Ba{x)Bp{x'). In particular, they are relevant for the asymptotic behavior of th e fun ctions 
{Ba{x)Bp{x')^{xi) ■ ■ ■ $(a;„)) for x x' . Of course, these exponents also appear in the representations ( 5.10 ) if the 
correlation function F{r) in the left hand side involves the operators with j > n. 



VIII. DISCUSSION AND CONCLUSIONS 



The zero-mode and RG techniques have been exploited in a model of magnetohydrodinamics turbulence where the 
magnetic field is passively advected by a Gaussian white- in-time velocity, in the presence of a const ant ba ckg round 
magnetic field that introduces a large-scale anisotropy. The basic equations of the model are Eqs. (pTl|)-( |2!^ ). We 
have shown that the correlation functions of the magnetic fluctuations exhibit inertial-range anomalous scaling. The 
explicit asymptotic expressions for the correlation functions of the magnetic field and their powers have been obtained. 
In the inertial range, the correlation functions are represented as superpositions of power laws with universal exponents 
and nonuniversal amplitudes. The anomalous exponents have been calculated both nonperturbatively (for the second- 
order correlation function) and perturbatively (for the second and higher-order correlation functions), in the first order 
of the exponent ^ and in any space dimension d. 

In the language of the zero-mode techniques, anomalous exponents are associated with scale invariant functions 
which are annihilated by the inertial operator L (remember that in defining £ we neglected the molecular diffusivity 
Ko): the so-called zero modes of the equations for the correlation functions. In the language of the RG, these exponents 
are d eter mined by the critical dimensions of tensor composite operators built of the magnetic field without derivatives, 
Eq. (6.4), and exhibit a kind of hierarchy related to the degree of anisotropy: the less is the rank, the less is the 
dimension and, consequently, the more important is the contribution to the inertial-range behavior. The leading terms 
of the even (odd) structure functions are given by the scalar (vector) operators. 

For the pair corr elatio n function, the complete set of the exponents has been calculated nonperturbatively using 
the exact equation (3.14); t hey are giv en in Eqs. (4.27) together with the discussion of their admissibility. 

The general expressions ( ^.10| ), ( |7.3|) describe the behavior of the correlation functions for Ar >> 1, and any fixed 
mr (m= 1/L) where A~^ (x r], ij being the dissipative scale, and L is the integral scale of the problem; expressions 
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( [7.2[ ), (7.4) correspond to the additional condition mr « 1 (inertial range). These results for the leading terms can 
be summarized as follows: 



Pn + q (X f 



with A„^p given by [see Eqs. ( |6.15 ) and ( |6.18 )] 

2p{p- 1) - (d- l)in-p){d + n + p) 



\x — X 



2(d- l)(d + 2) 



(8.1) 



(8.2) 



Here is some component of B, e.g., its projection onto the direction r/r or B°/B°, pn is the minimal possible 
value of p for a given n (i.e. p = for n even and p = \ for n odd). The exponents C"''' a-re expressed through the 



dimensions A„^p as follows: 



A 



n.O 



A,,o = - 



(d + 2) ^ 
^ {nq + rf 



o(C') 



if n, g are even, 



(d + 2) 



+ 0(^^) if n, q are odd. 



.3) 



(d + 2) 



Oie) 



if n is even and q is odd. 



In the presence of an anisotropic forcing, questions about isotropy restoration at small scales are naturally raised. 
In particular, an issue recently addressed concerns the behavior of the derivative skewness factor of the passive scalar 
at large Peclet number, Pe, in the presence of large-scale anisotropy, and, in a more general formulation, the effects 
of large-scale anisotropy on the inertial-range statistics of passively advected fields [^||^, [l^ , [l6|] and the velocity itself 
|p|, p^p^ . In the case of passive advection of a scalar field, both the real Q and the numerical experiments HJl^ show 
that the derivative skewness remains 0(1) for very high Peclet, in disagreement with what could be expected on the 
basis of both dimensional argumentations and cascade ideas. It means that, contrary to K41 hypothesis, anisotropy 
present at large scales persists at small scales. For the velocity field, in the case of an homogeneous shear flow an 
guivalent result has been found for the vorticity, which keeps a constant value, independent of the Reynolds number 

Let us now briefly discuss the consequences of our results for anisotropic indicators in this problem. Since the 
equation (^) is not invariant with respect to the shift B ^ B + const, we can use as the simplest measure of 
small-scale anisotropy the dimensionless ratios of the correlation functions of the field B without derivatives, e.g.. 



R„ 



{B--\x)B\\{x')) 
[B.{x)Bn{x'))-/^- 



From Eqs. (S^) it then follows that in inertial range of scales we have: 

R2k+i oc (Ar)-'^^'^'.o (TOr)'^^'=+i'i-(2fc+i)'^2-o/2^ 



(8.4) 



(8.5a) 



R2k+2 oc (Ar)-^^'=+i'i (mr) 



A2fc + 2,0-(fc + l)A2.0 



.5b) 



Note that the ratios ( |8.5D depend on both scales of wavenumber A and m; the dependence on the former follows 
from the fact that the powers have nontrivial anomalous dimensions. The dependence on the Peclet number, 

Pe = (A/m)^ can be estimated by replacing r with rj = 1/A; see Ref. ||l^. Using explicit O(^) expressions for A„_p 
we then obtain: 



(8.6a) 



i?2fe+2 0cPe2fe(fc+l)/(d+2). 



.6b) 



Since the leading terms of the even functions ( ^.l| ) are determined by the expon ents of t he is otropic shell (i.e., those 
related to scalar RG operators), the inertial-range behavior of the even ratios ( ^.5b| ), ( ^.6b| ) is the same as in the 
isotropic model. This gives a quantitative support to the unive rs ality o f anomalous exponents with respect to different 
classes of forcing. On the other hand, the odd quantities (B.5a), (S.6a) appear to be sensitive to the anisotropy: i?3 in 
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( B.6a ) slowly decreases for Pe oo, while ratios i?2fc+i with k > 2 increase with Pe. Moreover, general expressions 
( B.5a ) contain large A dependent factors, which also prevent these functions from vanishing at Pe oo. Notice the 
important difference between the isotropic and the anisotropic problem: in the former the (nonuniversal) constant of 
the inertial range power laws of odd-order moments are zero by symmetry, while this is not the case in the anisotropic 
context. This implies that the (hierarchical) exponents for the odd-order moments appear solely in the anisotropic 
case. For a given odd order, the leading exponent is thus responsible of the observed scale-dependent normalized odd 
order ratios. 

The picture outlined above seems rather general. Indeed it is compatible with that recently estab lished for the NS 
turbulence and for the scalar field passively advected either by the velocity of the type ( |2.2D (see Ref. |Q) or 

by a NS velocity in the two-dimensional inverse cascade regime (see Ref. ) . For a passive scalar field advected by a 
rapidly changing velocity field such as (2.2), RG expressions for the dimensionless ratios i?„ = Sn/S!^^^, Sn being the 
n-th order structure function of the scalar field, are given by the expressions ( ^.5| ) without A dependent factors and 
with the same exponents of mr (see Ref. for 6*3 and Ref. for general n's and d's). So, for example, for k = 1 
the ratio S2k+i decreases down to the small scales, but much slower than it was expected on the basis of dimensional 
argumentations, while for fc > 1 it grows in agreement with the results of Ref. 

It should be emphasized, however, that the results obtained within the lowest-order approximations in ^ are reliable 
only for moderate n, because the actual expansion parameter in the Kraichnan model is rather than ^ itself (see 
Ref. [p5| ). The analysis of the large n behavior requires some additional resummation of the ^ series, which remains 
an open problem. For the passive scalar case, the numerical experiments and the instanton calculus show 



that the exponents analogous to C"'^ in Eq. (ST) tend to a finite limit as n ^ 00 ("saturation"). It is worth noting 



that the limiting expressions for n ^ 00 obtained in Refs. | ]36| diverge as ^ — > 0, thus signalling that small ^ and 
large n limits do not commute. The persistence of small-scale anisotropics and the intermittency saturation are both 
statistical signatures of quasi-discontinuities observed in the scalar field . It is then reasonable to expect that the 
saturation of intermittency takes place also in the magnetic case where quasi-discontinuous structures in the magnetic 
field are likely to be present (see, e.g., |37y3^ ]). 
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APPENDIX A: COEFFICIENTS IN THE EQUATIONS FOR THE MAGNETIC CORRELATION 

FUNCTIONS 



We report hereafter the coefficients a^, bi, ■ ■ ■, ri appearing in Eqs. (4.6)-(4.9) 





1. Coefficients in 


Eq. (4.6) 


ai 






bi 


^id-i){d-i-0 
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TOi = -^(2C^ - 8C + 8 - 2d) 
m = 2^ 

01 = - 2) 

Pi = -C(C-2)(C-4) 

91 = - 2) 



2. Coefficients in Eq. ( ^ ) 

a2 = (d + e-l)[2-2^ + ^(^-l)] 
62 - d - 1 

C2 = (d+e-i)(rf-i) + 2e 
^2 = d + e - 1 

e2 = -(rf + e-l)(rf-l) 

/2 = + e(2d - 3) + _ 3d] 

.g2-2eK' + e(d-2)-d] 

m2 = -((i + C- 1) 
n2 = -(d + ?-l)(?-2) 



3. Coefficients in Eq. ( [4.8| ) 

as = (2-0(^^ + ^-1) 
fe3 = -e(rf + e-2) 
C3 = d - 1 

d3 = e + (d-i)' 
/3--(d + e-i)(rf+i) 

33 -d^ + 2d - .^d + 4^ - 1 - 2C^ 

J3 = -d^ 

fc3=e 

?3 = e(rf+e-2) 
"13 = -'2m - 2) 

= 2m - 2) 



4. Coefficients in Eq. ( |4.g| ) 

a4 = 2(2-0(rf + e-l) 
64 = d - 1 

C4 = (d + e- l)(d - 1) + 2^ 
d4 = d + ^ - 1 

e4 = -(d + C-l)(c'-l + 20 
/4 = -2^ 
54 = 2^ 
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APPENDIX B: RELATIONS INVOLVING THE LEGENDRE POLYNOMIALS 

From the well-known relations involving the Legendre polynomials Pj{z) (see, e.g., Ref. |^^) the following decom- 
positions for a function F{r, z) = X^jlo Pji'^)fji^) hold: 
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APPENDIX C: FULL SET OF EQUATIONS PROJECTED ON THE LEGENDRE POLYNOMIALS 



Inserting Eqs. and ( 4.13| ) into Eqs. ( |4.6| )-(4.11) and exploiting the relations reported in Appendix the 

following equations follow from the orthogonality of Legendre polynomials: 
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TABLE I. Canonical dimensions of the fields and parameters in the model (3.1) 
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FIG. 1. Behavior of zero-mode e xpone nts Q (j — 0, 2,4 and 6) vs ^ for d = 3. Notice the inequality Co < C2 < C4 ' 
means the validity of the hierarchy (4.22) and thus the self-consistency of our zero-mode solutions. 
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